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Abstract 

We determine to order a 3 the so-called residual mass in the lattice regularisation 
of the Heavy Quark Effective Theory for Nf = 2. Our (gauge-invariant) strategy 
makes use of Numerical Stochastic Perturbation Theory to compute the static 
interquark potential where the above mentioned mass term appears as an additive 
contribution. We discuss how the new coefficient we compute in the expansion of 
the residual mass can improve the determination of the (MS) mass of the 6-quark 
from lattice simulations of the Heavy Quark Effective Theory. 



1 Introduction 



Given our present computer resources, one can not accommodate the 6-quark on the lattice 
keeping at the same time under control both finite size and finite lattice spacing effects. 
Despite this fact, the determination of the 6-quark mass is one of the most precise results 
in Lattice Gauge Theory. There are various strategies that have been successfully applied 
to this purpose. A first determination comes from Non-Relativistic QCD [1]. A second one 
relies on Heavy Quark Effective Theory (HQET) [2] . Both these approaches have to deal with 
the mass counterterm whose computation we will be concerned with in this paper. In order 
to avoid this task - which has to be performed in Perturbation Theory (PT) - a third strategy 
has been devised, which makes use of HQET, but in a Non-Perturbative framework. A direct 
computation of the mass counterterm is avoided by making use of a finite scaling technique 
[3]. A fourth method has been recently introduced which actually avoids HQET and remains 
within QCD. It again relies on a finite scaling technique, from which it takes the name of 
Step Scaling Method [4]. A different approach - which can also profit of our computation 
- has been followed by [5]. When different strategies are available, an important issue is 
the comparison of their results. In the quenched approximation, all the previous methods 
are consistent within errors. At the moment the only unquenched computation is the one 
in [2]. However the largest source of error, there, comes from the limited knowledge of the 
perturbative residual mass. This can be strongly improved by the a 3 computation which we 
are presenting in this paper. One important issue concerning the mass counterterm, whose 
perturbative computation we are going to present, is that it is a power-divergent quantity. 
Because of that, there is no continuum limit in PT, and that was the motivation for the 
approach in [3]. Of course, while this is a conceptually important point, it is not a no-go 
for the first two methods. One has to stick to a finite value of the lattice spacing and to 
carefully assess the error which is inherent to this procedure. All these considerations set 
the stage for our computation. Our main task is to make it possible to extend the result in 
[2] to a higher accuracy in the unquenched case. The quenched determination has already 
been greatly improved by a previous quenched computation of ours [6] of which the present 
work is an extension. Having reached the goal of a higher accuracy in the analysis of [2], 
one can also investigate the status of PT with respect to the fact that one has to stick to a 
finite lattice spacing. 

The structure of the paper is as follows. In the second section we sketch the theoretical 
framework, while in the third we present our computational strategy. More details about 
those issues in common with the quenched computation can be found in [6]. The fourth 
section collects our results together with a brief discussion of the impact of the computation. 



2 The residual mass in Lattice HQET 

The Lagrangian density of HQET [7] may be written as 



where D± is the (gauge) covariant time derivative and h is the bi-spinor describing the only 
degrees of freedom of the heavy quark that are relevant in this approximation. While the 



£hqet — hD 4 h, 




1 



theory is appealing also in force of its simple form, it is well known that there are many 
subtleties one has to face, many of which have to do with a proper definition of the heavy 
quark mass itself. This is well evident in the fundamental relation one would like to exploit in 
order to connect the mass of a physical hadron (M B ) to the HQET expansion mass parameter 
(mj) and the (linearly divergent) binding energy (£) 



M B = m b + £ + 0(l/m b ). 



(2) 



At this stage m b is actually not yet properly defined. Following [2] (and references therein) 
one can proceed as follows. By matching the QCD propagator to its lattice HQET counter- 
part one gets a relation involving the pole mass 



m 



pole 



M B -£ + Sm + 0(l/m b ). 



(3) 



In this relation a new character has entered the stage: a linearly divergent additive mass 
counterterm, the so-called residual mass Sm which in lattice regularisation - contrary to 
dimensional regularisation - is generated by quantum corrections. Our task is now to compute 
this quantity in Perturbation Theory 



+i 
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n>0 



The pole mass (by writing m v h ole one has in mind the pole mass for the 6-quark) can be 



related to the MS mass m b 
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The perturbative matching reads 



m b (m b ) = m p b 



pole 
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(5) 



The coefficients D n with n < 2 are known [8, 9]. One can now put everything together to 
get what the authors of [2] call the master equation 



m b (m b ) 



M B -S + J2(a s (mE)) 
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(6) 



After having arrived to this fundamental relation a few comments are in order. Note, first 
of all, that the last expression asks for both a non-perturbative (£) and a perturbative 
(Sm) computation. As it is well known, £ can be computed as the decay constant of the 
correlation function of two axial currents. In Eq. (6) many things are actually taking place, 
for which the perturbative expansion of the residual mass is responsible. First, Sm is in 
charge of canceling the linear divergence of £. Also, a subtle cancellation of renormalon 
ambiguities is taking place: the renormalon in the expansion of Sm cancels the one in the 
perturbative relation between the pole and MS masses [10]. It should now be stressed that 
both things are taking place in Perturbation Theory. Because of that, the cancellation of 
the linear divergence is incomplete and the continuum limit can not be taken. As for the 
renormalon ambiguities cancellation, the same coupling must be used in both perturbative 
expansions present in Eq. (6). As we have already said, D 2 is known, which means that X 2 
(from which X 2 is computed) has to be computed in order to exploit at best the knowledge 
coming from continuum Perturbation Theory. X has been known for a long time, while Xi 
was computed in [11]. As already said, X 2 was computed in the quenched approximation in 
[6] and then also in [12]. Our goal is now to compute X 2 for Nf = 2 Wilson fermions, in 
order to make contact with the unquenched data of [2]. 
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3 The computational strategy 



As already said, the computational strategy is just the extension to Nf = 2 of the computa- 
tion in [6]. Consider a gauge invariant loop such as a Wilson (W) or a Polyakov (P) loop. 
By following [13] 's notation, one can write 

(W) = exp(-cL/a)W log . (7) 

In the above formula L is the length of the loop and a the cutoff scale at small distance (of 
course one can think of the lattice spacing), c is the linear divergence we are interested in. 
At the time of [13] it was regarded as a mass renormalization of the test particle that one 
considers if one reads the loop as an effective amplitude. The notation W\ og reminds that 
only logarithmic divergences are left. A first one is the divergence that can be absorbed in 
the definition of the renormalized coupling. A second logarithmic divergence appears if the 
contour has an angle; this is usually referred to as the corner divergence. Of course such a 
problem is there in the case of a Wilson loop and not for a Wilson line. One can now go to 
Lattice Perturbation Theory (LPT) and compute the static potential via Creutz's ratios 

V(R) ee Mm V T (R) V T (R) ee log (^^) ■ ( 8 ) 

The corner divergences now disappear because of the ratio, the number of corners being 
the same in numerator and denominator. As for the coupling, there is a standard way to 
renormalize it: 

V(R) = 25m + V C oui{R) V Coul (R) = ~C F ^^-. (9) 

This equation amounts to the definition of the coupling in the potential scheme. Since we 
are computing in LPT, what we are actually computing in the second formula of Eq. (9) is 
the matching between the potential and the lattice couplings. This reads 

V(R) = 2 5m - ^ (a + c^R) a 2 , + c 2 (R) ag + . . .) (10) 
where the relevant coefficients are given by 



Cl (R) = 2b \ogR + 2b \og^ 



Av 



c 2 (R) = c 1 (R) 2 + 2b 1 \ogR + 2b 1 \og-^ + 2 - 2 . 

A b 

The matching coefficients have been written down in terms of the dependence on A param- 
eters and coefficients of the /3-function. The good message is now that both the matching 
between the potential and the MS couplings ([14]) and the matching between the MS and 
the lattice couplings ([15]) are known. This means that one can match the potential and the 
lattice couplings and as a result the only unknown quantity in Eq. (10) is 5m. This defines 
our strategy: we computed V(R) in LPT up to the third order and fitted (order by order) 
5m to make contact with Eq. (10). 
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Figure 1: Comparing V^ 3 ^/?) at iVj = (squares) and Nf 
order of the potential with the expansion in terms of atQ. 



2 (circles). This is the third 



It is worth to comment that the only difference with respect to the procedure in [6] comes 
from the dependence on Nf of ci(R) and C2{R), which in turns comes from the dependence 
on Nf of ^ and ~ ^2' '• The final effect of this can be seen 1 in Fig. (1). 



3.1 Unquenching the computation 

Our computational tool was Numerical Stochastic Perturbation Theory (NSPT). This a 
stochastic way of performing Lattice Perturbation Theory which works by expanding the 
solution of the Langevin equation as a formal power series in the coupling 

k=l 

Such a solution has to be plugged in the Langevin equation 

^-U = {-iVS[U]-i V )U, 

1 Both in this figure and in the following ones results are plotted in terms of an expansion in ao- In [6] 
the plots were in terms of 1//3. 
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r\ being a gaussian noise and V being the Lie derivative on the group that acting on the 
action produces the classical equation of motion. In the limit of the stochastic time t going 
to infinite, expectation values with respect to the gaussian (77) noise reconstruct the results 
of expectation values with respect to the path integral measure [16]. Needless to say, the 
only difference in computing Wilson loops in the quenched and the unquenched case are 
the background gauge configurations, which should be computed with a new weight in the 
functional integral, i.e. 

e Sa ^ e S G det M = e -(S G -Tr InAO 

Sq being the Wilson gauge action and M being the lattice Dirac operator. The unquenched 
version of NSPT requires in the Langevin equation the evaluation of 

VS G - VTr InM = VS G - Tr ((VM)ilf- 1 ). 

The second term of the former expression can now be evaluated by a new source of noise 

(Re^VM)^ 1 ^)^, 

in which a sum is understood on repeated indices. This is an old strategy for unquenched 
simulations [17], which in NSPT gets greatly simplified by a simple observation. In our 
method everything should be expanded as a power series in the coupling, which means that 
also the matrix M* -1 ) has to be computed as the inverse of 

fc>0 

Now the inverse of this expression reads 

M -l = M (0)-! _ /3 -V2 M (0)- 1 M (l) M (0)-l + 

The form of the first non trivial order is enough to explain where the good point is. The 
only inverse to be computed is the zeroth-order, which is diagonal in momentum space and 
independent of the gauge field. The unquenched version of NSPT was first discussed in [18]. 
The computation of Wilson loops on which this work relies is its first actual application. 
It is worth to inspect again Fig. (1). The Nf = data are obtained from a sample of 120 
configurations 2 , while those at Nf = 2 from a sample of 200. As one can see from the 
errorbars, the statistical fluctuations are practically insensitive to the presence of fermion 
loops. A reader interested in some extra details on the method can refer to [19]. The 
computation of Wilson loops is actually only one out of a variety of computations for which 
the configurations we used were generated for (other goals have to do with full QCD, and in 
particular mainly with renormalization constants). 

4 Results 

4.1 5m to third order for Nf = 2 

To perform the present computation we thermalized 200 configurations on a 32 4 lattice with 
Wilson gauge action and Nf = 2 Wilson fermions, which is the unquenching that is relevant 

2 Actually these configurations are split into two sets with different discretizations of the time step in the 
Langevin equation. This is needed in order to perform a continuum stochastic time extrapolation. 
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to supplement for an extension of the results in [2]. 100 configurations were produced at 
dt = 0.01 (time step discretization of the Langevin stochastic time), and 100 at dt = 0.02. 
The extrapolation errors at dt = are included in the errorbars in Fig. (1). The mass of 
both sea quark was set to zero, by exploiting the knowledge of k cr u, which - at this order - 
has been analytically computed in [20, 21]. We computed W(R,T) for all R, T < 16 up to 
order o^. Again, a reader interested in the actual results for Wilson loops can refer to [19]. 
Before presenting our result, we now discuss in some more detail our procedure, together 
with some comments on errors. 

• First of all, statistical errors are associated with the numerical nature of NSPT. These 
are magnified by the process of extracting the potential, which requires computing a 
series out of the direct expansions one obtains for the primary quantities (which are 
the Wilson loops). The statistical errors associated to our measures were computed 
via the bootstrap technique [22]. 

• Given the finite lattice size nature of our computation, we can not actually take the 
T — > oo limit in Eq. (8) and it is clear that Eq. (10) is distorted by lattice artifacts. 
Given the overall finite size and the presence of two scales (i.e. R and T/R) it follows 
that there is no simple finite scaling technique one can apply within the constraints 
of our computation. One issue is of course to take both R and T/R fairly large. One 
reassuring piece of information concerning lattice artifacts came from the fitting of the 
constants coefficients entering Eq. (10) (i.e. Cp, 2& log^ il , . . . ). If one lets them as 
free parameters, their range of variation within the chosen value of R and T/R (see 
next point) is only of a few percent. 

• We determined the preferred values of the Xi by fitting our expansions of Vr(R) to 
Eq. (10). In these fits the Xi were the only free parameters (i.e. all the other constants 
were assigned their infinite volume, continuum limit values). The intervals in which 
fits were performed were such that three conditions were fulfilled: R > 3, T > 12 
and T > 2.5R (we denote by R the mean value of R in the fitting interval; the fitting 
intervals themselves ranged from 3 up to 7 points). It was reassuring to find out that 
under these conditions good x 2 values were found. 

• As for the values of \ 2 m the fitting intervals we referred to in the previous point, a first 
observation is that they are quite high for the fit of Xo: the best values are of order 6, 7. 
This is not surprising, given our fixed lattice artifacts approach: measurements of 

are simply precise enough to clearly detect lattice artifacts effects. For the coefficients 
Xi and X 2 we quote an error that embraces the range of values for which \ 2 does not 
exceed the value 2 (again, in the same fitting intervals in which both R and T / R were 
fairly large). 

Fig. (2) and Fig. (3) show an example of our fits at each order. 

The final results which we obtained are X = 2.118(3) (to be compared with the analytical 
X^ = 2.118...), X~i Nf=2) = 10.58(4) (to compare with [11] X^ {Nf=2) = 10.588...) and finally 

X~ 2 (Nf=2) = 76.7(6). (11) 
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Figure 2: A typical result of the fitting procedure described in the text. The plots refer to 
first and second order. The solid lines are the expressions of Eq. (10) once 6m has been 
fitted, while circles are the Vt{R)- In this case T = 14. Expansions are in terms of oq. 



4.2 The impact of the result 

The impact of the quenched result of [6] has been known since the Lattice 2000 conference 
[23]. As for the unquenched case, the authors of [2] have now repeated their analysis for the 
6-quark mass taking into account the new result for X { 2 Nf ~ 2) . This results in [24, 25] 

m^(m&) (Unque) = (4.21 ± 0.03 ± 0.05 ± 0.04) GeV 

The last error is the one taking into account the indeterminations in the perturbative match- 
ing, on which our result has an effect. It would be twice as large with no knowledge of X 2 
(with a central value for the mass of 4.26; see [23]). The new analysis [24] pins down different 
values for the quenched and unquenched results, even if they are compatible within errors, so 
that further investigation is needed. As for the finite lattice spacing dependence (remember 
that in this approach one sticks to finite values of the lattice spacing), it is not dramatic and 
it gets decreased by including the new term in the matching. On the other hand, the control 
on renormalon ambiguities seems quite firm. 
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Figure 3: The same as Fig. (2) for the third order. Again, the expansion is in terms of o.q. 



5 Conclusions 



We computed the perturbative expansion of the unquenched (Nf = 2) residual mass term in 
lattice Heavy Quark Effective Theory to order ajj. This allows to improve the determination 
of the 6-quark mass from (unquenched) lattice simulations of the HQET. Errors appear 
to be under a fairly good control and the impact on the final result for the quark mass is 
quite important, even if further investigation is needed if one wants to carefully assess the 
difference between quenched and unquenched results for the 6-quark mass. 



Acknowledgments 



We are very grateful to V. Gimenez which on behalf of all the authors of [2] has shared with 
us the results of their analysis after our new result has been taken into account. Our interest 
in the subject was first triggered some years ago by G. Martinelli and C.T. Sachrajda, from 
both of which we learned much on the subject. F.D.R. acknowledges support from both 
Italian MURST under contract 2001021158 and from I.N.F.N. under i.s. Mill. L.S. is 



8 



supported by DFG through the Sonderforschungsbereich 'Computational Particle Physics' 
(SFB/TR 9). 



References 

[1] See for example S. Collins, "The mass of the b quark from lattice NRQCD," in 
Quark Confinement and the Hadron Spectrum, World Scientific, 325 (2002). arXiv:hep- 
lat/0009040. 

[2] V. Gimenez, L. Giusti, G. Martinelli and F. Rapuano, JHEP 0003 (2000) 018 
[arXiv:hep-lat/0002007]. 

[3] J. Heitger and R. Sommer [ALPHA Collaboration], JHEP 0402 (2004) 022 [arXiv:hep- 
lat/0310035]. 

[4] G. M. de Divitiis, M. Guagnelli, F. Palombi, R. Petronzio and N. Tantalo Nucl. Phys. 
B 675, 309 (2003) [arXiv:hep-lat/0305018]. 

[5] G. S. Bali and A. Pineda, Phys. Rev. D 69 (2004) 094001 [arXiv:hep-ph/0310130]. 

[6] F. Di Renzo and L. Scorzato, JHEP 0102, 020 (2001) [arXiv:hep-lat/0012011]. 

[7] M. Neubert, Phys. Rept. 245, 259 (1994) [arXiv:hep-ph/9306320]. 

[8] K. G. Chetyrkin and M. Steinhauser, Phys. Rev. Lett. 83, 4001 (1999) [arXiv:hep- 
ph/9907509]. K. G. Chetyrkin and A. Retey, Nucl. Phys. B 583, 3 (2000) [arXiv:hep- 
ph/9910332]. 

[9] K. Melnikov and T. v. Ritbergen, Phys. Lett. B 482, 99 (2000) [arXiv:hep-ph/9912391]. 

[10] G. Martinelli and C. T. Sachrajda, Phys. Lett. B 354 (1995) 423 [arXiv:hep- 
ph/9502352]. 

[11] G. Martinelli and C. T. Sachrajda, Nucl. Phys. B 559, 429 (1999) [arXiv:hep- 
lat/9812001]. 

[12] H. D. Trottier, N. H. Shakespeare, G. P. Lepage and P. B. Mackenzie, Phys. Rev. D 65, 
094502 (2002) [arXiv:hep-lat/01 11028]. 

[13] V. S. Dotsenko and S. N. Vergeles, Nucl. Phys. B 169, 527 (1980). 

[14] Y. Schroder, Phys. Lett. B 447, 321 (1999) [arXiv:hep-ph/9812205]. 

[15] C. Christou, A. Feo, H. Panagopoulos and E. Vicari, Nucl. Phys. B 525, 387 (1998) 
[Erratum-ibid. B 608, 479 (2001)] [arXiv:hep-lat/9801007]. 

[16] G. Parisi and Y. s. Wu, Sci. Sin. 24, 483 (1981). 

[17] G. G. Batrouni, G. R. Katz, A. S. Kronfeld, G. P. Lepage, B. Svetitsky and K. G. Wilson, 
Phys. Rev. D 32, 2736 (1985). 



9 



[18] F. Di Renzo and L. Scorzato, Nucl. Phys. Proc. Suppl. 94, 567 (2001) [arXiv:hep- 
lat/0010064]. 

[19] F. Di Renzo and L. Scorzato, JHEP 0411 (2004) 036 [arXiv:hep-lat/0408015]. 

[20] E. Follana and H. Panagopoulos, Phys. Rev. D 63 (2001) 017501 [arXiv:hep- 
lat/0006001]. 

[21] S. Caracciolo, A. Pelissetto and A. Rago, Phys. Rev. D 64 (2001) 094506 [arXiv:hep- 
lat/0106013]. 

[22] R. Alfieri, F. Di Renzo, E. Onofri and L. Scorzato, Nucl. Phys. B 578 (2000) 383 
[arXiv:hep-lat/0002018]. 

[23] V. Lubicz, Nucl. Phys. Proc. Suppl. 94, 116 (2001) [arXiv:hep-lat/0012003]. 

[24] V. Gimenez, private communication. 

[25] P. E. L. Rakow, arXiv:hep-lat/0411036. 



10 



